Checks of asymptotia in pp elastic scattering at LHC 



Agnes Grau a , Simone Pacetti b , Giulia Pancheri c , Yogendra N. Srivastava b 

" Departamento de Fisica Teorica y del Cosmos, Universidad de Granada, 18071 

Granada, Spain 
b INFN & Physics Department, University of Perugia, Italy 
^vq C INFN Frascati National Laboratories, Via E. Fermi 40, Italy 

o 

(N 

Abstract 

We parametrize TOTEM data for the elastic differential pp cross section at 
y/s = 7 TeV in terms of two exponentials with a relative phase. We employ 
& two previously derived sum rules for pp elastic scattering amplitude in impact 

^ parameter space to check whether asymptotia has been reached at the LHC. 

^ A detailed study of the TOTEM data for the elastic differential cross section 

at \/s = 7 TeV is made and it is shown that, within errors, the asymptotic 
sum rules are satisfied at LHC. We propose to use this parametrization to 
i—i study forthcoming higher energy data. 
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1. Introduction 

O 

CN| The TOTEM experiment has measured the total pp cross section at 

^/s = 7 TeV [1] and the value may be consistent [2] with saturation of 
the Froissart-Martin bound J3j 0], i.e. cr total ~ In 2 s. It is also consistent 

^ with a range of total cross section predictions, based on a mini-jet QCD 

model with infrared gluons [51 E] , as well other deductions based on a hard 
Pomeron component [7j. The saturation of the Froissart bound, if confirmed, 
could exclude, at least at present energies [8j, a different behaviour obtained 
from hidden extra dimensions [91 [10]. Indeed, according to Ref. [11], the 
TOTEM result indicates that we have now reached asymptotia as far as the 
total cross section is considered and thus existence of hidden extra dimen- 
sions can be excluded. On the other hand, contradictions in the analysis 
of Ultra High Energy Cosmic Rays composition [12] might hint to the ex- 
istence of new physics that modifies baryonic interactions at center-of-mass 
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(CM) energies in the 50 -j- 100 TeV range. Thus the problem of the energy 
dependence of the total cross section is not yet beyond all dispute. 

Another interesting feature has arisen in connection with TOTEM pub- 
lished results for the differential and the total elastic cross section [13]. These 
data show that the dip, which had characterized data at ISR in pp scatter- 
ing, and which has been only a faint presence in pp at the TeVatron, has 
reappeared. We propose to use these data to check whether asymptotia has 
been reached, through two asymptotic sum rules (ASR) for the elastic am- 
plitude in impact parameter space [H]. To check the satisfaction of these 
ASR, we shall use the model independent parametrization of the scattering 
amplitude [U] proposed in 1973 by Barger and Phillips. 

We shall see that this model, with two exponentials and a phase, describes 
well both the diffraction peak, as well as the t-dependence after the dip at 
LHC with CM energy of 7 TeV (LHC7). We also analyze ISR data for pp 
scattering at y/s = 53 GeV, and, by comparison, find that the ASR are not 
satisfied at ISR, while being almost exactly satisfied at LHC7. We discuss 
the approach to asymptotia of the ratio cr e i as tic/ctotai an d of the forward slope 
and suggest to use this simple model to fit future forward scattering LHC 
data. 

2. Asymptotic sum rules for the elastic scattering amplitude 

Let the elastic amplitude F(s, t) be normalized so that 



At high energy, all particle masses can be ignored and the elastic differential 
cross section reads 



ctotaits) = 47r F(s, 0) . 



(1) 



da 



n\F(s,t)\ 2 . 



(2) 



Writing the elastic amplitude as 



Jo 

inversion of Eq. (|3j) reads 



/■oo 

F{s,t) = i / (bdb)J (bV^i) [l - e ^«(M e -2MM] ; 



(3) 
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and one can study properties of the amplitude in impact parameter space 
from the known behaviour in the angular momentum discrete variable. The 
Froissart bound implies that there must exist a finite angular momentum 
value, below which all partial waves are absorbed. Under the stronger 
hypothesis of total absorption as b — > in the ultra high energy limit, namely 
all "low-6" waves are completely absorbed, we have the two ASR 

1 f° 

SRi = - / (dt)^smF(s,t) — > 1; as s — > oo , (5) 

^ J-oc 

i r° 

SR = - (dt)9teF(s,t) ->■ 0; as s^ oo . (6) 

^ J-oo 

Satisfaction of these rules is a good measure of whether the asymptotic limit 
has been reached, and would reinforce the statement [2] based on the TOTEM 
data for total cross section, that we may have reached asymptotia, namely 
that crtotai ~ In 2 s, and saturation of the Froissart bound has been observed. 
Notice that the Froissart-Martin bound is obtained under a hypothesis weaker 
than total absorption [I] and it would lead to SRi —> 2; SR —> 0. As we 
shall see, phenomenologically Eq.([5]), is favoured. 

In order to check these ASR, one requires a model for the scattering am- 
plitude. While the imaginary part of the elastic amplitude is solidly anchored 
to the optical theorem, the real part is more model dependent. A limiting 
behaviour for the real part of the forward amplitude has been obtained by 
Khuri and Kinoshita [161, i.e., for s — > oo 



MeF(s,0) 7T 
P(S ' } ~ QmF(s,0) ~ Ms/ s ) ' [ ) 

if the Froissart bound is saturated. As discussed later, the model of Ref. [TT] 
gives p(s, 0) = 7r/[2pln(s/s )] with p = 1 for a total cross section asymptoti- 
cally growing as ~ In s, and p = 1/2, when the Froissart bound is saturated. 

For values of q 2 = —t ^ 0, there is a result by Andre Martin [I8l [19], 
which relates the real part of the amplitude to its imaginary part. Given an 
5smF(s,t), Martin's method gives an asymptotic expression for $leF(s,t), 
i.e. 

3teF(s,t) = p(s,Q)^-[tQmF(s,t)]. (8) 

We note here that the above equation is valid also if the Froissart bound 
is not saturated, but o" to tai ~ ln 1 ^ p (s), where 1/2 < p < 1. The elastic 



3 



amplitude still scales with the variable r = ta tota i, the only difference is 
p = 7r/[2pln(s/s )], with s = 1 GeV 2 . 

To complete the expression for the real part, one needs to estimate p(s, 0) 
for all cases including the case when the Froissart bound is not saturated. 
The phase of the leading (first term) contribution at t = is readily obtained 
from the s-dependence of the total cross section, or equivalently QmF(s,0), 
using the prescription s — > se~ ln ^ 2 , which assures analyticity and crossing 
symmetry of the scattering amplitude [T8l 120] . At high energies and very 
small t values, s -H-tt crossing symmetry tells us that the leading C = + 
amplitude is a function of the complex variable se~ ,7r / 2 . The Froissart bound 
requires QmF(s, 0) < In 2 s, and we can generally write, asymptotically 



0) -> 
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with l/2<p<lto include both the case of saturation of the Froissart 
bound and a slower rise, compatible with Ins behaviour. For large s, this 
may be approximated to 
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2p\n(s/s )_ 



(10) 



The above equation gives an estimate for the leading contribution to the 
parameter p(s,0), namely 

3?eF(s,0) 



p(s,0) 



7T 



(11) 



3mF(s,0) 2pln(s/s ) ' 

where the approximate sign refers to the fact that the real part of the am- 
plitude can receive also (non-leading) contributions from other terms. The 
Khuri-Kinoshita result, valid when the Froissart bound is saturated, is ob- 
tained for p = 1/2. To obtain a numerical value for p, we use a model j6] 
built on QCD mini-jets embedded into the eikonal representation. Accord- 
ing to this model, QCD mini-jets drive the rise of the cross section, and k t - 
resummation of ultra-soft gluons transforms the power-law rise of the mini-jet 
cross sections into an asymptotic logarithmic behaviour. In our model [TT] 

a totaJ ~ ln^{s) , (12) 

where the parameter 1/2 < p < 1 can be related to a dressed one-gluon 
exchange potential rising as r 2p_1 , and our phenomenology of total cross 
section data [5] leads to values p = 0.66 -r 0.77. 
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It should be noted that the asymptotic expression for the real part of the 
scattering amplitude automatically satisfies the second ASR. However, this 
does not help in checking whether asymptotia has been reached, since the 
expression for the real part is based on an asymptotic scaling law, which may 
not be valid. 

Another measure for asymptotia, also examined in Ref. [2], can be ob- 
tained from the ratio of the elastic to the total cross section. Based on the 
black and grey disk model, where Ototai = 2tt AR 2 (s) and elastic = 7iA 2 R 2 (s), 
the black disk limit, A = 1, gives, for the ratio, the limiting value 1/2. How- 
ever, we have recently shown [B] that the eikonal formulation at high energy 
underestimates the total inelastic cross section, including in it only uncorre- 
lated (non-diffractive) collisions, whereas the diffractive collisions appear to 
be part of the eikonal elastic cross section. In fact, it is known that, at high 
energies, the "elastic" cross section, as defined in eikonal models, is typically 
overestimated ETj. Thus we propose the following asymptotic limit 



sr-, ^"elastic ~r ^diffractive t> , t> ,1/0 , ZioN 

7c = = 7c e i + TCdiff — V 1/2 , s — > 00 , (13) 

°"total 

where (Jdiffractive includes both single and double diffractive components. We 
show in Fig. [T]a compilation of data [22] for the ratio cr e i ast i C /cT tc .tai and com- 
pare it with the black disk limit and results of Ref. (2). In the figure we 
have added to the accelerator data a value extracted from the recent mea- 
surement by the Auger Collaboration for the inelastic cross section [23J. 
We estimated the ratio lZ e \ at 57 TeV, by using Block and Halzen (BH) [2] 
value for the total cross section at cr™ al (57 TeV) = (134.8 ± 1.5) mb, which 
is based on the analytic amplitude method of Ref. [21]. We then obtained 
a c iastic(57 TeV) = cr t B H al - = (44.8 ± 11.6) mb. We also show the 

asymptotic result (green dot) from Ref. [2], which appears in line with the 
ratio (Xeiastic/cTtotai extrapolated from the Auger point for the inelastic cross 
section. The conclusion in Ref. [2J that we are still not in asymptotia for 
what concerns this ratio, may however be too pessimistic. Firstly, from the 
figure, it is quite conceivable that data may approach a limiting value of 



1/3. This would be better in keeping with Eq. (13). Also, estimates for the 
inelastic cross section are affected by large uncertainties, and often model 
dependent on extrapolations. Indeed, in order to reconcile cosmic ray com- 
position data, an extrapolation from p-air to pp could use higher total cross 
section values 1121. Given the uncertainties in different simulations at ultra 



high energies, the estimate for TZ e \ at 57 TeV could then be lowered (closer 
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Figure 1: The ratio Ceiastic/ctotai compared with accelerator data (black triangles) |22j . 
including TOTEM. The value at y/s — 57 TeV is obtained from the ratio of the Auger 
collaboration value for cri ne i [23] to the value for ototai predicted in Ref. [2]. 



to the TOTEM value), thereby opening the possibility that an asymptotic 
plateau for creiastk/ctotai may have already been reached. In fact lZ e \ has been 
rising from lower energies until the value ~ 1/4 at the TeVatron, and at 
LHC7. Data at 8 TeV could indicate whether we have reached a plateau or 
not. 



3. The two exponential model at LHC energies 

In order to perform a test of the ASR in a model independent way, we 
have analyzed the data for the differential elastic cross section through a 
parametrization of the scattering amplitude, proposed by Barger and Phillips 
(BP) in 1973 [15]. 

In 1973, BP proposed two different parametrizations, consisting in using 
a phase and two exponentials to describe the elastic scattering amplitude as 
a function of momentum transfer in the range —t = 0.15 -f- 5.0 GeV 2 . The 
first parametrization was given as 

A(s, t) = % \y/A(8je? B W + y/CWeWe&W] , (14) 
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^ = A(s)e B ^ t + C{s)e D W + 2y/A{s)y/C{s)e^ +D ^ 2 cos <\> , (15) 

with five s-dependent real parameters, A, B, C, D, (p. The total cross 
section in this model is given as 

ototai = 4v^r3wt(s,£ = 0) = \/A(s) + y/C(s) cos0 . (16) 



Using the above expression, BP fitted data from pi ab = 12 GeV to 1496 GeV 
(y/s w 5 -T- 53 GeV). The values of the parameters at yfs = 53 GeV are 
given as VJ = 6.55, \[C = 0.034 in ^mb/GeV 2 and B = 10.20, D = 1.7 in 
GeV -2 , = 2.53 rad. The first exponential was seen to have normal Regge 
shrinking, namely 

B = B + 2a Ins, a' » 0.3 , (17) 

while the second exponential term appeared to be constant in this energy 
range. With this parametrization, the phase was fitted to be always 
larger than n/2, so that the interference term was always negative. The 
second parametrization interpreted the energy dependence of the parameters 
in terms of Regge-Pomeron exchange and fits to data were as good as in 



the first case, but more model dependent. However, Eq. (17) for the for- 



ward slope B poses now a problem, as recently discussed in Ref. [23], since 



TOTEM data can be described by Eq. (17) only if a' is not a constant, but 



has a logarithmic energy dependence. On the other hand, by relinquishing 



any model interpretation, one can use the model of Eq. (14) as it appears 
to be the simplest possible way to describe the t-dependence of pp elastic 
scattering amplitude from ISR onwards. 



The simplicity of the model in Eq. (14) suggests to use it to describe 
present LHC data. Reading the TOTEM data from Ref. [13], we show the 
result in the left panel of Fig. [2} We find that two exponentials and a phase 
can give a very good description of TOTEM data, with the tail well described 
by a term e Dt , with B/D ~ 4 and a relative weight of the two terms at LHC7 
yj A/C ~ 22. With this parametrization, the position of the dip (which 
the model fit gives at —t = 0.52 GeV 2 ) and the behaviour around the dip 
are both well reproduced. This result differs from the proposal in Ref. [7] 
and also from the suggestion by TOTEM [13] that the tail, past the dip at 
—t = (0.53 ± 0.01 stat ± 0.01 syst ) GeV 2 , is described as being compatible with 
a behaviour of the type \t\~ n with n = 7.8 ± 0.3 stat ± 0.1 syst . The TOTEM 
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suggested behaviour is close to a description of the tail through a QCD type 
contribution from three gluon exchange advanced in Ref. [26J. It is also 
close to the observation in Ref. [15], that the "second exponential" can be 
identified with a term proportional to p^ M , valid, according to the authors, 
for all available pp data for s > 15 GeV 2 . Following these suggestions, and 
mindful of the fact that the form factor dependence of the amplitude would 
contribute to the cross section with a term as (— 1)~ 8 , we have tried a slight 
modification of the BP model, which consists in substituting the second term 
in Eq. (14) with a term proportional to the 2 nd power of the proton form 



factor, namely to write 



A(s,t) = i 



VW)e Bt/2 + 



V 7 ^ 



(-t + t y 



(18) 



With to a free parameter, this allows a comparison with the parametrization 
of the tail in Ref. [13], but the fit worsens, and the tail is everywhere better 
described by an exponential. We show this fit in the right hand panel of 
Fig. [2} From here on, we shall use the two exponential model and proceed 



s=(7 TeV)" 



A = 458 ± 5 mb GeV" 
B = 17.21 ±0.06 GeV" 2 
C = 0.90 ± 0.03 mb GeV" 2 
D = 4.38 ± 0.03 GeV" 2 
<|> = 2.70 ± 0.01 rad 
Z 2 /d.o.f. = 2.25 
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Figure 2: Five-parameter fit of pp data for the elastic differential cross section at LHC7. 
At left the fit uses the two exponential model from Ref. |15j . at right a form factor term 



is used to describe the tail, according to Eq. (18). Parameter values in the two cases are 



indicated. TOTEM data were read from Ref. 113 . 



to check the ASR. Since the phase 7^ tt, for the satisfaction of both ASR 
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the model of Eq. ( 14 ) needs to be modified so as to include a real part in the 



leading amplitude iy/A(s), i.e. the leading term (at t ~ 0) should acquire 
an extra phase. Rewriting the amplitude of Eq. (fl4j) as 



one has 



%Ai + A 



R 



e \B{s)t + 



i \/C(s) cos <p — VC sin i 



p(s,0) 



A h 



C sin< 



Aj + \[C cos ( 



(19) 



(20) 



We use the result for the leading contribution to p(s, t — 0) to give a real 
part to the first term, i.e. 



p(s,t = 0) 



A f 



71 



Aj 2p]n(s/s y 



(21) 



With this, we now find 



SRq 



A 



C sin i 



1 + p 2 



B' 



D 



0.045 0.066 



(22) 



for p = 0.77 and 0.66 respectively, and 



SRi 




COS ( 



0.94 



(23) 



with A = A 2 R + A 2 -, and inserting the parameter values from the figure, with 
7r/2 < <p < 7r. One observes a good satisfaction of both ASR. Notice however, 



that, in general, the conjecture for p(s, 0) proposed in Eq. ( 11 ) cannot be valid 
at lower energies, because the approximation for the amplitude in Eq. ^ is 
only valid at asymptotic energies. Thus, for uniformity with the ISR case, 
curves in Fig. |2j have been done with Ar = 0. Applying the above model 
to the ISR data for pp scattering at y/s = 53 GeV [27], we obtain the fit 
shown in the left panel of Fig. [3] We notice that the application of the ASR 
to pp scattering at these lower energies, gives SRi ~ 0.7 for the parameter 
values obtained by the fit. As expected, the ASR for the imaginary part of 
the amplitude in 6-space is not yet satisfied at ISR. As for the second ASR at 
ISR, the expression for the real part of the amplitude as, discussed in Sec. 2, 
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Figure 3: At left, the description of pp scattering data at ISR, using the model of Eq. ( 15 ) 



A five-parameter fit gave the indicated values for the parameters. At right we show the 



i-dependence of the effective slope at LHC7 from Eq. ( 26 ) 



cannot be applied here, because at ISR we are far from an asymptotic regime. 
Before leaving this discussion, we point out that the value obtained for the 
total cross section in this model is below the experimental value. This is 
particularly true for ISR energies, where the fit gives cr to tai = (35.6±1.7) mb to 
be compared with values with very small errors, such as (42.71±0.35) mb [28J. 
At LHC7 the total cross section is better reproduced by the fit, with a value 
of (91 ± 5) mb, although still a bit low. From Figs. |2]and[3j we conclude that 
this model is approximately valid at lower energies, better apt to describe 
data at LHC7 and onwards. 

4. The slope parameter in the two exponential model 

TOTEM released values for crtotai, ^elastic, the slope B exp and the position 
of the dip. The total cross section is reproduced by the five-parameter fit 
(^totai is let unconstrained) of the two exponential model within 10%, and 
we have a good determination of the position of the dip, which, in our fit, 
occurs at — tdi P = 0.52 GeV 2 . 

As for the slope parameter, one cannot immediately compare the value 
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of this model parameter B with TOTEM data, where it is defined as 



dt dt 



t=o 



and is measured in two t-intervals, with B exp = (20.1 ±0.2 stat ±0.3 syst ) GeV 2 
for smaller \t\ values, i.e. 0.02 < -t < 0.33 GeV 2 , and B exp = (23.6±0.5 stat ± 
0.4 s y st ) GeV" 2 from the interval 0.36 < -t < 0.47 GeV 2 . Such difference is 
well understood in the two exponential model. With the definition for the 
effective slope [20] 

B«{8,t) = ±h(^) , (25) 



dt \ dt 

one obtains a function changing with t, i.e. 

„ , , ABe Bt + CDe Dt + y/AVC(B + D)e( B+D W 2 cos <j) 

= da-Jdt ' (26) 

Using the parameter values from the left hand panel of Fig. [2j we show in 
the right-hand panel of Fig. [3] the t-dependence of the slope B e s at LHC7, 
compared with the values given by the TOTEM experiment in the two 
different t intervals, whose extension is indicated by the horizontal bars. At 
t — 0, one finds 



BoS (7TeV,0 H A ln (dg! 



= 17.7 GeV" , (five-parameter fit) .(27) 

t=o 



Notice that the figure is obtained by neglecting the contribution from Ar. 
Clearly this approximation is not valid where the real or the imaginary part 
of the amplitude approaches zero. It is possible to identify the two t-values 
where this will happen before the dip, namely where the real part of the 
amplitude is zero, and a subsequent one, close or very close to the dip, where 
the imaginary part is zero. At LHC7, these two points occur at 

= 0.28 + 0.30 GeV 2 , (28) 




2 ■ ( ^AW) \ 



In f ' = 0.5 GeV 2 , (29) 



B — D \ |cos0| I 
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where the uncertainty in i R depends on the value for p and we have neglected 
p 2 in Eq. (29) and, again, we have used the parameter values indicated in 



the left hand panel of Fig. [2j 

An asymptotic value for the slope -B e ff(s,i), can be obtained through the 
ASR, already well satisfied at LHC7. Thus, in the exact limits SRq = 0, 
SRi = 1, for asymptotic energies at LHC and beyond, one obtains for the 
effective slopes 

B cS (s,t R ) = ^, (30) 
B^{s,t I ) = B + D- a ^. (31) 

We find, in agreement with TOTEM, that past tn, where the real part of the 
amplitude is zero, the slope is ~ 20 GeV -2 , and that it increases to higher 
values, to then start decreasing as it approaches tj where the imaginary part 
goes to zero in correspondence with the occurrence of the dip. 

We note that the ASR imply that the slope B e g(s,tji) from the two ex- 
ponential model should grow as er to tai, thus as In 2 s, if the Froissart bound is 
saturated, or as ln 1/,p s for a slower rise. We see that the same behaviour is 
also true asymptotically for the (constant in t) leading parameter B of the 
two exponential model, for which we also have the well known result 



/ \ ""elastic ) V -A "total /QO^ 

recently also discussed in Ref. |29j on general grounds. 

We conclude that, if the two exponential model is a good representation 
of the amplitude at LHC in the range < -t < 3 GeV 2 , and if the total 
cross section rises faster than In s, the slope of the leading term in this model 
should grow faster than m(s/so). Notice that in Ref. [22], an analysis of the 
effective slope from NA8 to TOTEM seems to indicate a rise according to 
In 2 s. Thus, just as in the case of the total cross section, the question of the 
energy dependence of the effective slope is still very much open. 

In conclusion, we stress the following points concerning the slope param- 
eter: 

1. One can not speak about a unique slope if there is more than one term 
with different t dependence in the elastic amplitude, as required by the 
presence of a dip, in the differential cross section. 
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2. Even in the region before the dip, the necessity of at least two terms, 
whose interference gives rise to the dip, is felt and TOTEM itself pro- 
vides two different values for the slope in two different regions of t. 
Hence it is not meaningful to consider only a single slope and its asymp- 
totic behavior. In a two component model such as the one considered 
above, all one can discuss is the relative growth of the two slopes with 
energy. 

3. Even in the rather simple model analyzed here, we find that the slope 
near t = is less than that around t = —0.2 GeV 2 , which is less than 
that around t = —0.4 GeV 2 and hence the question of the growth with 
energy becomes a strong function of the t range under consideration. 

The curve in the right hand panel of Fig. [3] shows that at LHC this model 
predicts an increase of the slope parameter in the range < — t < 0.4 GeV 2 . 
This is in contrast with the behaviour observed at lower energies, notably at 
ISR for very small t— values [301 E] • Notice however that our findings at LHC 
are in agreement with those expected in Ref. |32j at asymptotic energies. 

Conclusion 

We have utilized two asymptotic sum rules derived earlier for the elas- 
tic scattering amplitude in pp scattering [H] and have tested them using a 
well known parametrization of the elastic differential cross section with two 
exponentials with a relative phase. Our study of TOTEM data at LHC7 
finds that this model can describe well the data before and past the observed 
dip. Similar conclusions about the behaviour past the dip, have also been 
advanced in Ref. |33j . 

A possible interpretation of this model is that the exponential behaviour 
exhibited both before and after the dip corresponds to a resummation of soft 
terms accompanying the leading C = +1, two gluon exchange, and C — — 1, 
three gluon exchange term, with non leading terms giving rise to <fi ^ it. 

Given the general nature of the parametrization discussed here, we 
recommend its usage for further studies of the behaviour of the scattering 
amplitude at higher energies. 
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